A new limit of T-violating transverse muon polarization in the Tx^ji^v decay 
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A search for T-violating transverse muon polarization (Pt) in the K^ ■k^^'^v decay was 
performed using kaon decays at rest. A new improved value, Pt ~ —0.0017 ± 0.0023{stat) ± 
0.0011{syst), was obtained giving an upper limit, \Pt\ < 0.0050. The T-violation parameter was 
determined to be Im^ = -0.0053 ± 0.0071(stat) ± 0.0036(si/st) giving an upper limit, |Im^| < 0.016. 



PACS numbers: 11.30.Er, 12.60.Fr, 13.20.Eb 

The transverse muon polarization, Pt, in the 
TT^/j,"*"!^ (^/ts) decay is one of the observables of CP vio- 
lation beyond the standard model (SM) . CP violation in 
general is a subject of continuing interest in K and B me- 
son decays. Pt, defined as the polarization component 
perpendicular to the decay plane, is an obvious signature 
of a violation of time reversal (T) invariance, since the 
spurious effect from final state interactions is very small 
(< 10"^) Q. Pt is almost vanishing 10"'^) in the SM 
with the Kobayashi-Maskawa scheme |^ ; it is therefore a 
very sensitive probe of CP violation mechanisms beyond 
the SM and new physics along with B physics such as 
6 ^ S7 and some other decays |3J- Models 3 such as 
those with multi-Higgs doublets or leptoquarks, or some 
SUSY may to give rise to Pt as large as lO^^^. 

At the High Energy Accelerator Research Organiza- 
tion (KEK) the E246 collaboration has been performing 
a search for Pt in K'^^. In 1999, the first result was pub- 



lished ^ based on ^ 3.9 x 10^ good K^^^ events from the 
data taken during 1996 and 1997, indicating no evidence 
for T violation. Further runs provided a cumulative data 
sample with three times more events. This Letter consti- 
tutes our final result from all the data with an improved 
analysis. The present result supersedes all our earlier 
reports. 

The principle of the experiment was the same as de- 
scribed before 0. A kaon beam with an average inten- 
sity of 1.0 X 10^/s was produced at the 12 GeV pro- 
ton synchrotron from 3 x 10^^ protons per spill of 0.7 
s duration with a 2.7 s repetition time. The detec- 
tor setup (Fig.l) using stopped kaons at a 12-sector 
magnet is described in detail in 6]. The muon po- 
larization consists of three components a) longitudinal. 
Pi — 'TV/IPmI parallel to the muon momentum p^, 
b) normal, Pn ^ s^- {p^ x (pV x p^))/|p^ x (pV x p^)| 
normal to p^ in the decay plane, and c) transverse. 
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FIG. 1: Experimental setup, (a) cross section side view, (b) end view, and (c) cross section view of the muon polarimeter of 
one sector at a certain radial position r with the tilted positron counters. The y direction is also shown. See for details. 



Pt = S/i • {Ptt X Ph)/\'Ptt X P^L\ perpendicular to the de- 
cay plane. Pt was searched for as the azimuthal polar- 
ization (0 or y component in Fig.l) of /i^ emitted radi- 
ally (in the r direction) and stopped in the Al stoppers 
when a tt" was tagged in the forward (fwd) or the back- 
ward {bwd) direction relative to the detector axis. The 
spin depolarization during flight and in the stopper was 
estimated to be negligible. This azimuthal polarization 
was measured as an asymmetry between clockwise (cw) 
and counter-clockwise (ccw) emitted Michel , New and 
Nccw Summation over the twelve sectors with 12-fold az- 
imuthal symmetry played an important role in reducing 
systematic errors. Events from fwd and bwd tt'^s have 
opposite asymmetries. We exploit this feature to double 
the signal and also as a powerful means to cancel the 
systematic errors. 

The total data were grouped into three periods of (I) 
1996-1997, (II) 1998, and (III) 1999-2000, each having 
nearly the same beam conditions and amount of data. As 
was described in j5||, two completely independent analy- 
ses, Al and A2, pursued their own best off-line event se- 
lections with their own analysis policies. This approach 
provided a cross-check of the quality of selected events 
and also an estimate of the systematic errors in the anal- 
ysis. Basic event selection criteria were same in both 
analyses. The 7r°'s were identified not only as two pho- 
tons (27) but also as one photon (I7) with energy >70 
MeV. The maximum sensitivity to Pt is provided by the 
fwd and bwd regions of 7r° (27) or photons (I7) with 
I cos 6'jrO(-y) I > 0.342, where 6'jrO(-y) is the polar angle cor- 
rected for muon direction. Slight differences between the 
two analyses led to a non-negligible amount of uncom- 
mon good events in each analysis. All the selected events 
were categorized into the common {Al ■ A2) events and 
two sets of uncommon events (^1 • A2 and Al ■ A2) sepa- 
rately for 27 and I7. In total 6.3 million and 5.5 million 
good events were obtained for 27 and I7, respectively. 
The fraction of 27 and I7 mismatch events between Al 
and A2 was only 1.5% and these were rejected. The 



positron yield was extracted from the time spectra by 
integrating from 20 ns to 6.0 fis after subtraction of the 
constant background deduced from fitting between 6.0 
fis to 19.5 /is. The only significant background to muon 
stopping and its decay was due to tt"*" decay in flight from 
if+ 7r+7r°; its contamination effect was estimated and 
included in the systematic errors. 

In the T-violating asymmetry At was calculated 
as At = {Rf/Rb - l)/4, where = {Ncw/Nccw) f{b) 

for the TT^-fwd {bwd) region, using the total positron 
cw and ccw counts. Then, Pt was calculated as Pt = 
AT/{aint < cos 9t >) using an average analyzing power 
aint and the angular attenuation factor < cos 6t > with 
9t being the angle of decay plane normal vector relative 
to the y axis. However, this method is prone to a sys- 
tematic error due to potentially different muon stopping 
distributions of fwd and bwd events. To obtain a finite 
stopping efficiency, muon stoppers with finite size in the y 
and r directions were employed. A geometrical asymme- 
try appears for muons at {y, r) off center which, in turn, 
can induce a fake At if the muon stopping distribution 
is different between fwd and bwd events, in particular in 
the y direction. In the current analysis an exact treat- 
ment, in which we use the y muon stopping point from 
the C4 tracking chamber located just in front of the stop- 
per, was employed. For the r direction an integration 
was used because the change of geometrical asymmetry 
is much smaller (about 1/10 of y-dependence), and be- 
cause its determination from tracking was poor. The 
transverse polarization Pt for each data set was evalu- 
ated as the average of contribution PT{y) from each part 
of the stopper using the C4 y coordinate as; 

Pt= I PT{y)w{y)dy (1) 



where w{y) is the weight function proportional to 
l/(Tp^(y) (here, ap^{y) is the error distribution) and nor- 
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malized to 1, and Pxiv) is 



At{v) 



a{y) < cos9t > 



(2) 



with the y-dependent asymmetry Axiy) and analyzing 
power a(y). The definition of At (y) = [{Af{y)—Ai,{y)]/2 
assured that it was free from the geometrical asymmetry 
and from muon stopping densities, and canceled the sys- 
tematic errors common for fwd/bwd. Here, Af(y) and 
Ab{y) were calculated as Af(^t){y) = [{Ncw{y)/Nccw{y) - 
l)/2]^^£;(-b^^). The y dependence of analyzing power 
could be calibrated using the positron asymmetry AN{y) 
associated with the normal polarization as AN{y) is 
proportional to a(i/). Apf was measured by rearranging 
the fwd and bwd events into left and right categories of 
7r° directions, and calculating An — [Aieft — Aright) f'^^- 
This has a maximum at the center of the stopper (Fig. 2). 
The absolute value of a was calibrated by a Monte Carlo 
simulation as a = Aff'-^/P^'-'. The coefficient a{y) in- 
cluded the effects of intrinsic muon decay asymmetry, 
muon spin precession around the field, positron interac- 
tions, and the finite counter solid angle. The obtained 
a{y) function corresponded to ~ 0.271 ± 0.027, 

which is significantly higher than our previous estimate 
of a^nt = 0.197 ± 0.005 ,5] deduced as a = A^^^^/P^'^ 
and thus less rehable. Pt thus obtained in Eq.(l) is re- 
garded as the average value of Pt distribution in the 
finite kinematical acceptance of if^a in the stopper. The 
validity of applying the proportionality relation Priy) ^ 
AT{y)/AN{y) in Eq.(2) was carefully checked under the 
actual trigger condition. In order to increase the sta- 
tistical accuracy of a{y), An of all the data sets was 
summed since a(?/) is only dependent on y and should 
not depend significantly on the data set. In the actual 
analysis, the averaging of +i and —i bins was used be- 
cause the shape of «(?/) should be symmetric in the first 
order approximation also in the presence of the magnetic 
field. Fig. 2 shows Priv) thus calculated which is nearly 
constant with slight but opposite-sign gradients for 27 
and I7. This is due to different muon stopping distribu- 
tions along the r direction for fwd and bwd events with 
an opposite tendency of S < r >=< r >fwd — < r >bwd 
for 27 and I7 due to kimematics. Pt was calculated, for 
the integration Eq.(l), by summation over the 36 bins 
from y = —9.0 cm to -1-9.0 cm. The effect of the r-origin 
Priy) gradients is eliminated since the effect cancels be- 
tween +y and —y, and the y distribution is symmet- 
ric. The average values of y, weighted by the statistical 
significance of respective PT{y), were < y >=0. 007mm 
and 0.020 mm for 27 and I7, respectively. These small 
< J/ >'s confirm an excellent C4/stopper alignment and 
justify this analysis. The factor < cos^t > was evalu- 
ated for each data set by using a Monte Carlo calculation 
taking into account realistic background conditions, to be 
typically 0.7 and 0.6 for 27 and I7, respectively. 

Data quality checks were performed for the 18 data 
sets of the 3 groups with 6 data categories each. First 
the null asymmetry was calculated as the asymmetry of 
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FIG. 2: Measured An (upper) and Pt (lower) as functions of 
y. Black dots (•) are 27 events and open circles (o) are I7 
events. Each bin is 0.5 cm wide with the center (y = cm) 
at between 18 and 19 bins. 



all the fwd and bwd events added, using the total cw 
and ccw counts integrated over y and it was confirmed 
that there was no significant bias. Next, An were com- 
pared. Although there was a slight difference among the 
I7 data sets due to different cut criteria of the event 
selection we decided to use all the I7 data. Then the 
distribution of decay plane normal {n^o x ) with its 
Or and Oz components was studied to check for any 
possible kinematical phase space distortions, and no sig- 
nificant offsets were found. Finally, the 18 Pt values 
(Table I) which are consistent with each other (a fit to 
a constant gives /v = 0.78, where v is the degree of 
freedom), yielded the average of Pt = -0.0017 ± 0.0023 
being consistent with zero. The sector dependence of 
Pt is plotted in Fig. 3 with /v = 0.69 for 27 data 
and x^/v — 1.97 for I7 data showing that the latter 
is slightly inferior. The Pt's were converted to the T- 
violating physics parameter Im^ (Tj] with the conversion 
coefficients $ = 0.327 and 0.287 for 27 and I7, respec- 
tively, deduced from a Monte Carlo simulation The 
ideogram of Im^ (Fig. 4) shows that there is significant 
overlap among the different data sets. The average is 
found to be Im^ = -0.0053 ± 0.0071. It is noteworthy 
that the analysis by the previous method gives consistent 
central values of Pt = -0.0018 and Im^ = -0.0063. 

Although almost all the systematics were canceled due 
to the summation of the 12 sectors and the double ratio 
between fwd and bwd events, a few errors remain giving 
rise to spurious At or a small admixture of Pn resulting 
in a spurious Pt effect (Table II). The contribution of 
misalignments of detector elements and the muon spin 
rotation field remained as in 'B^ . The small mean values 
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TABLE I: T-violating polarization Pt of the 18 data sets of 27 and I7 events from the two analyses of Al and A2 for three 
experimental periods of I, II and III. The errors are only statistical. For the definitions of event categories, see the text. 



Data category 



1(1996-1997) 



11(1998) 



111(1999-2000) 



27 [Al ■ A2] 
27 [Al- A2] 
27 [Al •A2] 
I7 [Al ■ A2] 
I7 [Al- A2] 
I7 [Al -^"2] 



0.00112 ±0.00667 
-0.00735 ± 0.01022 
-0.00385 ± 0.00899 
-0.01393 ± 0.00956 
0.01014 ±0.01069 
0.00228 ±0.01134 



-0.00317 ±0.00729 
0.01225 ± 0.00858 
0.00640 ± 0.01268 
-0.01366 ±0.01042 
-0.01114 ±0.01280 
-0.01660 ±0.01531 



-0.00596 ±0.00711 
-0.00037 ± 0.00754 
-0.00473 ± 0.01201 

0.01113 ±0.01035 
-0.01088 ±0.01022 

0.00951 ±0.01195 
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FIG. 3: Dependence of Pt on the sector number. Black dots 
(•) are 27 events and open circles (o) are I7 events. 
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FIG. 4: Ideogram of Im(5). Black dots (•) are data sets I, 
open circles (o) are II, and stars (★) are III. /v = 0.78. 



of 9r and 9z were treated as an error. The effect of muon 
multiple scattering through the Cu degrader may cause 
a difference in the actual muon stopping distribution of 
fwd and bwd, in particular in the y distribution even for 
a measured y at C4, producing a spurious At through 
the geometrical asymmetry along y. This effect, inad- 
vertently omitted in our previous analysis, was carefully 
estimated in the present analysis to be SPt — 7.1 x 10^'*. 
The small effect due to Priy) gradients and finite < y > 
values was treated as a systematic error {SPt = 2 x 10~^) 
and included in the item "Analysis" together with other 
analysis uncertainties . The total systematic error was 
calculated as the quadratic sum of all the contributions 
resulting in APt — 1.1 x 10"'^ which is much smaller 
than the statistical error. 

In conclusion, we obtained the values of Pt — 



-0.0017 ± 0.0023(stat) ± 0.0011{syst) and Im$ = 
-0.0053 ± 0.0071 (stat) ± 0.0036{syst) with no indication 
of T violation. The 90% confidence limits are given as 
\Pt\ < 0.0050 and |Im$| < 0.016 by adding statistical and 
systematic errors quadratically. This result is a factor 3 
improvement over the last BNL-AGS experiment 8] and 
it may constrain the lightest Higgs mass and/or other 
parameters in the framework of non-SM models 4] bet- 
ter than or complementary to the neutron electric dipole 
moment c?„ and B meson decays. For example, our re- 
sult gives a stronger constraint to the three Higgs doublet 
model than the similar semileptonic decay B — > XtVt Q 
and implies in one of the multi-Higgs models (Q Garisto 
and Kane) that the down quark contribution to dn should 
be more than a factor 10 less than the current experimen- 
tal limit of d^^P < 6.3 X 10-2^ecm; our |Im^| 90% limit 
corresponds to 0.5 x 10~^^ecm. 
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TABLE II: Summary of systematic errors. 

Source SPt x 

e+ counter misalignment 2^9 

Misalignments of other counters 2.6 

Misalignment of B field on the stopper 6.1 

K'^ stopping distribution < 3.0 

Decay plane rotations {6r and 9z) 1.4 

/i^ multiple scattering 7.1 

Backgrounds (including tt^ decay from ) < 2.0 

Analysis (including Pt gradients) 4.0 



Total 



< 11.4 
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